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A complete, rigorous relativistic field-theory formulation of the nucleon-nucleon (NN) 
bremsstrahlung reaction is presented. The resulting amplitude is unitary as a matter of course 
and it is gauge invariant, i.e., it satisfies a generalized Ward-Takahashi identity. The novel feature 
of this approach is the consistent microscopic implementation of local gauge invariance across all 
interaction mechanisms of the hadronic systems, thus serving as a constraint for all subprocesses. 
The formalism is quite readily adapted to approximations and thus can be applied even in cases 
where the microscopic dynamical structure of the underlying interacting hadronic systems is either 
not known in detail or too complex to be treated in detail. We point out how the interaction cur- 
rents resulting from the photon being attached to nucleon-nucleon-meson vertices can be treated 
by phenomenological four-point contact currents that preserve gauge invariance. In an advance 
application of the present formalism [K. Nakayama and H. Haberzettl, Phys. Rev. C 80, 051001 
(2009)], such interaction currents were found to contribute significantly in explaining experimental 
data. In addition, we provide a scheme that permits — through an introduction of phenomenologi- 
cal five-point contact currents — the approximate treatment of current contributions resulting from 
pieces of the NN interaction that cannot be incorporated exactly. In each case, the approximation 
procedure ensures gauge invariance of the entire bremsstrahlung amplitude. We also discuss the 
necessary modifications when taking into account baryonic states other than the nucleon TV; in de- 
tail, we consider the A(1232) resonance by incorporating the couplings of the AW to the NA and 
AA systems, and the 7iV — s> A transitions. We apply the formalism to the 280-MeV bremsstrahlung 
data from TRIUMF [Phys. Rev. D41, 2689 (1990)] incorporating 7iVA transition currents and find 
good agreement with the data. 

PACS numbers: 25.10.+S, 25.40.-h, 25.20.Lj, 13.75.Cs 
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I. INTRODUCTION 

The two-nucleon system is one of the simplest strongly 
interacting systems. The study of the nucleon-nucleon 
(NN) bremsstrahlung reaction, therefore, offers one of 
the most fundamental and direct avenues for understand- 
ing how the electromagnetic field interacts with strongly 
interacting hadronic systems. In the past, the NN 
bremsstrahlung reaction had been applied extensively 
mainly to learn about off-shell properties of the NN in- 
teraction. It should be clear, however, that off-shell ef- 
fects are model-dependent and cannot be measured, and 
therefore are meaningless quantities for the purpose of 
comparison. 

Even though the original motivation for investigating 
the NN bremsstrahlung reaction has fallen away, under- 
standing the dynamics of such a fundamental process, 
nevertheless, is of great importance from a general theo- 
retical perspective. This is highlighted by the fact that 
none of the past models of NN bremsstrahlung could de- 
scribe the high-precision proton-proton bremsstrahlung 
data from KVI p], HJ for coplanar geometries involving 
small proton scattering angles. This was generally con- 
sidered all the more surprising since the irrelevance of 
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off-shell effects was taken as implicit proof positive that 
the coupling of a photon to the interacting two-nucleon 
system was under control. The discrepancy between the 
KVI data and the existing theoretical models, therefore, 
was quite unexpected. This longstanding discrepancy of 
nearly a decade was resolved recently by the present au- 
thors [3| who put forward a novel approach to the NN 
bremsstrahlung reaction that takes into account details 
of the photon coupling to interacting systems that had 
previously been neglected. The study, in particular, re- 
vealed the importance of accounting for the correspond- 
ing interaction currents in a manner consistent with the 
gauge- invariance constraint. 

Another recent bremsstrahlung experiment concerns 
the hard bremsstrahlung process p + p — > pp( So) + 7 
measured for the first time by the COSY-ANKE Col- 
laboration 0. In the absence of free systems of bound 
diprotons, this process was considered as an alternative 
to the 7+pp( 1 «S'o) — > P+P process which complements the 
photo-disintegration of the deuteron. Here, the hardness 
of the bremsstrahlung is due to the fact that the invariant 
mass of the two protons in the final state is constrained 
experimentally to be less than 3 MeV above its minimum 
value of twice the proton mass. In this kinematic regime, 
the two protons in the final state are practically confined 
to the 1 5o state and most of the available energy is car- 
ried by the bremsstrahlung. Therefore, this kinematic 
regime is as far away from the soft-photon limit as pos- 
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sible. The proton-proton (pp) hard bremsstrahlung re- 
action has been also measured at CELSIUS-Uppsala |5(. 
In spite of extensive studies of the NN bremsstrahlung 
reaction in the past, no dedicated experiments of pp hard 
bremsstrahlung with the diproton in the final state had 
been available until these recent measurements [!, 
Also, apart from the very recent study of Ref. Q, the 
theoretical investigation of the p + p — > pp^So) +7 reac- 
tion is virtually inexistent so far. 

Apart from the intrinsic interest in the elementary NN 
bremsstrahlung process, the investigation of this reaction 
has also an immediate impact in the area of heavy-ion 
physics. Indeed, di-lepton production in heavy-ion col- 
lisions is used intensively as a probe of hadron dynam- 
ics in the nuclear medium. Due to their weak interac- 
tion with hadrons, di-leptons are well suited to probe 
the hadron dynamics in the dense region of heavy-ion 
collisions. The HADES Collaboration, in particular, is 
currently engaged actively in the study of di-electron pro- 
duction in heavy-ion as well as in elementary NN colli- 
sions in the 1-2 GeV/u energy domain Q. In order to 
interpret the experimental data in heavy-ion collisions, 
it is imperative to understand the underlying basic el- 
ementary processes. Unfortunately, these basic elemen- 
tary processes are not yet fully under control. In fact, 
there is a number of theoretical efforts to understand 
these basic reaction processes Q . According to these 
studies, among the various competing mechanisms, the 
NN bremsstrahlung is one of the major mechanisms in 
producing di-electrons in these reactions. 

On the theoretical side, the majority of the exist- 
ing models of NN bremsstrahlung are potential mod- 
els. They have been applied to the analyses of exper- 
imental data (mostly in coplanar geometries) obtained 
up until the early 2000s, before the more recent experi- 
ments mentioned above were performed. Among those 
models, the most recent and sophisticated ones that 
have been used in the analysis of the high-precision KVI 
data [l], Bare the microscopic meson-exchange models 
of Refs. |10i4l2| . There are also a number of other mi- 
croscopic model calculations throughout the 1990s fl3l — 
l2~i| . which are dynamically similar to Refs. flB4l2j ]. ad- 
dressing a variety of issues in the NN bremsstrahlung 
process. All these models satisfy current conservation 
(at least in the soft-photon approximation) Q but none 
of them obey the more general gauge-invariance condi- 
tion in terms of the generalized Ward-Takahashi identity 
(WTI) employed in Ref. [|| (and explained in more detail 



The current conservation of earlier models usually comes about 
because for pp bremsstrahlung, there are no exchange currents 
for (uncharged) mesons and the four-point contact current dis- 
cussed in Ref. 3j is absent for phenomenological meson-nucleon- 
nucleon form factors that depend only on the momentum of the 
exchanged meson. For pn bremsstrahlung (see, e.g., Ref. [13|0 
the meson-exchange currents are taken into account via Siegert's 
theorem which preserves current conservation in the soft-photon 
limit. 



in the present work) . Quite recently, it was shown for- 
mally [22] how to maintain gauge invariance in a theory 
of undressed non-relativistic nucleons if one introduces a 
finite cutoff in a reference theory that is presumed to be 
already gauge invariant. When applied to effective field 
theories, in particular, this implies that gauge invariance 
can be maintained in such theories order-by-order in the 
expansion. 

The purpose of the present paper is to present the 
complete, rigorous covariant formulation of the NN 
bremsstrahlung reaction whose successful advance appli- 
cation was reported in Ref. The approach is based 
on a relativistic field theory in which the photon is cou- 
pled in all possible ways to the underlying two-nucleon T- 
matrix obtained from the corresponding covariant Bcthc- 
Salpeter-type NN scattering equation. This formulation 
follows the basic procedures of the field-theoretical ap- 
proach of Haberzettl [23| developed for pion photopro- 
duction off the nucleon. The resulting bremsstrahlung 
amplitude satisfies unitarity and gauge invariance as a 
matter of course. The latter, in particular, is shown ex- 
plicitly by deriving the corresponding generalized WTI. 

The formalism is quite readily adapted to approxima- 
tions and thus can be applied even in cases where the 
microscopic dynamical structure of the underlying inter- 
acting hadronic systems is either not known in detail or 
too complex to be treated in detail. As a case in point, we 
mention that the success of the advance application [3j of 
the present formalism to the KVI data [l|, [2| was due to 
the incorporation of phenomenological four-point contact 
currents that preserve gauge invariance following the ap- 
proach of Haberzettl, Nakayama, and Krewald [24[ based 
on the original ideas of Refs. [23|,[25j]. In addition, we pro- 
vide a scheme that permits the approximate treatment 
of current contributions resulting from pieces of the JVJV 
interaction that cannot be incorporated exactly. In each 
case, the approximation procedure ensures gauge invari- 
ance of the entire bremsstrahlung amplitude. We also 
discuss the necessary modifications when taking into ac- 
count baryonic states other than the nucleon iV; in detail, 
we consider the A(1232) resonance by incorporating the 
couplings of the NN to the NA and AA systems, and 
the 7iV — »■ A transitions. 

In Sec. HH we present the details of the full four- 
dimensional relativistic formulation, including a proof of 
the gauge invariance of the resulting bremsstrahlung am- 
plitude. We also introduce the necessary modifications 
for a covariant three-dimensional reduction and discuss 
its implications for the description of the dynamics of the 
process. We point out, in particular, that if one aims for 
a dynamically consistent microscopic description of all 
reaction mechanisms, one must implement gauge invari- 
ance in terms of generalized Ward-Takahashi identities 
for each subprocess — mere global current conservation 
is not sufficient. We show how one can preserve the gauge 
invariance of the amplitude even if some interaction- 
current mechanisms — both for hadronic three- and four- 
point functions — cannot be incorporated exactly. Fur- 
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thermore, we discuss what needs to be done to add ad- 
ditional baryonic degrees of freedom, in particular, the 
coupling of NN, NA, and AA channels. In Sec. HTJ 
we report on our results for the TRIUMF data [Hj]. A 
summarizing assessment of the present work is given in 
Sec. IIVI Some technical details follow in the Appendix. 

We emphasize that the present formalism is completely 
general and applies to proton-proton, proton-neutron as 
well as neutron-neutron bremsstrahlung processes. Fur- 
thermore, the photon can be either real or virtual. The 
former corresponds to the usual AW bremsstrahlung pro- 
cess while the latter case is suited for applications in di- 
lepton production in both the elementary NN and heavy 
ion collisions. 



II. FORMALISM 

The bremsstrahlung current B^ L is obtained from the 
nucleon-nucleon T-matrix by attaching an outgoing pho- 
ton to all reaction mechanisms of T in all possible ways. 
To this end, we use here the gauge-derivative procedure 
developed in Ref. [23[ in the context of pion photopro- 
duction. This procedure is formally equivalent to em- 
ploying minimal substitution for the connected part of 
the hadronic Green's function, and then taking the func- 
tional derivative with respect to the electromagnetic four- 
potential A 1 *, in the limit of vanishing A^ (for details, see 
[23| ) . The current is then obtained by removing the prop- 
agators of all external hadron legs from this derivative 
in a Lehmann- Szy manzik-Zimmermann (LSZ) reduction 
procedure [27l l29j. 

A. Full four-dimensional formalism 

The nucleon-nucleon T-matrix is determined by the 
corresponding four-dimensional Bethe-Salpeter scatter- 
ing equations, 

T = V+ VG T or T = V + TG V , (1) 

where V is the AW interaction given by the set of all two- 
nucleon irreducible scattering mechanisms. Go describes 
the intermediate propagation of two non-interacting nu- 
cleons, i.e., schematically we have 

Go = [h o ta] , (2) 

where U denotes the propagator of the individual nucleon 
i and "o" stands for the convolution of the intermediate 
loop integration. 

The bremsstrahlung current is obtained by evalu- 
ating the LSZ-type equation 

B^ = -G^ 1 {G a TG o y G^ 1 (3) 

where {• • • } M denotes the gauge derivative [23[ taken here 
of the connected hadronic NN Green's function GqTGq, 



G o d"G = 




FIG. 1. Graphical representation of Eq. (|8]l. Solid lines depict 
nucleons and wavy lines indicate the outgoing bremsstrahlung 
photon. 

with /i being the Lorentz index of the current. Using then 
the product rule {YXy = Y{Xy + {YyX for an (or- 
dered) product YX of a two-step sequence of hadronic re- 
action mechanisms described by operators X (first step) 
and Y (second step), we employ Eq. ([T]) repeatedly to 
find 

= (l + TGo)(d^ + V^)(G T+l)-d^ , (4) 

where d M defined by 

d" = -G^iGoVGo 1 (5) 

subsumes the one-body current contributions from the 
individual nucleons and 

v* 1 = -{vy (6) 

is the interaction current resulting from attaching the 
photon to any internal mechanisms of the NN interac- 
tion. Details of will be discussed below. 

Explicitly, the photon contributions from the two- 
nucleon propagator are found as 

d" = Tl(S 2 t^) + (Mr 1 )^ , (7) 

where is the electromagnetic current operator of nu- 
cleon i; Si denotes an implied S function that makes the 
incoming and outgoing momenta for the intermediate 
spectator nucleon i the same. We thus have 

G d"G = [hT^h o t 2 ] + [h o taT^ta] (8) 

which is represented graphically in Fig. [1] 

Note that Eq. (U) — apart from the subtraction by 
d M — possesses the structure of a distorted-wave Born 
approximation (DWBA), with the factors (GoT+ 1) and 
(1 + TGo) supplying the M0ller operators producing the 
initial and final scattering states, respectively, distorted 
by the NN interaction. The cP contribution by itself — 
without any initial-state or final-state NN interactions 

- is disconnected, as one sees clearly from Fig. [TJ The 
overall subtraction of d M in Eq. (UJ) , therefore, is necessary 
to remove this (unphysical) disconnected structure from 
B^ and retain only connected physical contributions. 

It is possible — and indeed desirable for the following 

— to equivalently rewrite Eq. Q to provide a full DWBA 
structure for i? M of the form 

B" = (l + TGo)>(G T+l) , (9) 

where 

J" = d>"GoV + VGod" + V 1 "- VG d^G V (10) 
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FIG. 2. Basic production current J M of Eq. (|lip for NN —¥ NN'y. Time proceeds from right to left. External legs are labeled 
by the four-momenta of the respective particles. Boxes labeled V subsume all two-nucleon irreducible contributions to the NN 
interaction that drive the Bethe-Salpeter equation (flj. The interaction current V* contains all mechanisms where the photon 
emerges from within the interaction V (i.e., any current mechanism not associated with an external leg of the interaction V). 
The last diagram, D M , as given in Eq. (|12p , subsumes all possible completely transverse contributions [cf. Eq. (|15|l ], in addition 
to the subtraction that corrects the double counting arising from the first two contributions (see text). Diagrams where the 
photon emerges from the lower nucleon line are suppressed; the antisymmetrization of nucleons is implied. 



is the completely connected current that describes the 
bremsstrahlung reaction in the absence of any hadronic 
initial-state or final-state interactions (with the exception 
of the subtraction VGod^GoV; see below); we shall re- 
fer to this tree-level-type current as the basic production 
current. The equivalence of this form of to Eq. (j4j is 
easily seen by repeated applications of Eq. ([1} . The NN 
T-matrices appearing to the right or left of the basic pro- 
duction current in © thus provide the initial-state 
interaction (ISI) or final-state interaction(FSI), respec- 
tively, of the two nucleons external to the basic produc- 
tion current J M . 

The subtraction term VGq(1^GqV in J' 1 removes here 
the double counting of contributions TGod^GoT to the 
full current B^ that arise from the two contributions 
in J M . As such, therefore, it is not a dynamically inde- 
pendent contribution to B^ and appears here only be- 
cause, for formal reasons, we wish to retain the DWBA 
form of B^ in Eq. ©. We shall see below, in Sec. Ill CI 
where we treat the covariant three-dimensional reduction 
of Eq. , that special considerations are needed for this 
subtraction term if one wants to maintain gauge invari- 
ance for the reduced amplitude. 

Equations © and (ITU1) are not the complete solution 
of the bremsstrahlung problem yet, because the gauge- 
derivative procedure — indeed any procedure based on 
minimal substitution — cannot produce current contri- 
butions that are completely transverse. Such contribu- 
tions must be added to the mechanisms obtained above 
for B^. Without lack of generality, we may do so by 
modifying the basic production current J M according to 



where 



J" = d"G V + VGod" + + D» (11) 



= T" - VGod^GoV (12) 



contains the sum of all explicitly transverse five-point 
currents denoted by T' 1 , in addition to the subtraction 
current VGad^GaV. In other words, 



Pi © — P2 



FIG. 3. Generic example for meson transition currents 
containing a transverse 7M1M2 vertex, with Mi 7^ M2, oc- 
curring in of Eq. (|12|l via its transverse contribution T M . 
The present numerical results were obtained incorporating 
7P7T and jujit transition currents. 



then is given by 



B" = (l+TGo)J»(GoT + l) 



(14) 



We emphasize that this equation and Eq. (jlip provide an 
exact generic description of the bremsstrahlung process 
off the NN system. The structure of the basic production 
current J M is depicted in Fig. [5J 

The detailed nature of the transverse contribution T 1 ' 
in (|12[) must be specified by the underlying interaction 



Lagrangians. Examples are meson-transition currents 
as depicted in Fig. [3] and 7iVA-transition currents 
J^, i.e., 



T n = j 



M 



(15) 







(13) 



The latter will be discussed below, in Sec. Ill FL when we 
consider A(1232) contributions in detail. 

The specific details of any particular application, of 
course, depend on the mechanisms taken into account 
in the NN interaction V that drives the scattering pro- 
cess in the Bethe-Salpeter equation (JTJ. For driving in- 
teractions based on single-meson exchanges, the com- 
plete structure of J M is discussed below, in Sec. Ill Dl 
In Sec. Ill Fl as mentioned already, we also consider the 
structures arising from A contributions that go beyond 
single-meson exchanges. 



B. Gauge invariance of the full four-dimensional 
amplitude 



is true irrespective of whether the external nucleons are 
on-shell or not. The complete bremsstrahlung current 



It should be clear that the procedure used for deriv- 
ing the bremsstrahlung current B^ does produce a cur- 
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rent that is gauge invariant as a matter of course. Nev- 
ertheless, we will now explicitly prove gauge invariancc 
because this will guide us later, in Sec. Ill C| in how to 
implement gauge invariance when we calculate in a 
covariant three-dimensional reduction. 

To prove the gauge invariance of the current (fT4|) . 
we first note that for an outgoing photon with four- 
momentum k, the four-divergence of gP may be schemat- 
ically written as 

M M = QGo 1 - G^Q , (16) 
where Q is short for 

= Qi + Q 2 . (17) 



The operator Qi describes the charge Qi of nucleon i 
and it removes the four-momentum k carried away by 
the photon from any subsequent interaction of nucleon 
i (appearing on the left of Qi). This notation allows 
one to keep track of the kinematics without explicit four- 
momentum arguments, i.e., specifying initial momenta 
for the two nucleons, the placement of the Qi immedi- 
ately allows one to find the momenta along each nucleon 
line. Equation (|16[) is an immediate consequence of the 
Ward-Takahashi identity for the nucleon current opera- 
tor Tf [HHl, i.e., 

W = Oiti 1 ~ t^Qi , (18) 

applied to Eq. ([7]). Explicitly, with Go specified as in ©, 
we have 



QGo 1 - g 'Q = KMr 1 - t^Qi] ° *2 1 + tl 1 ° - t^Oa] 

= [Qi^Cpi) - tiHpi - *0Qi] ° *2 l ip2) + t^ipx) o [Q 2 t 2 1 ( P2 ) - q 1 ( P2 - k)Q 2 ] (19) 

for a two-nucleon system where p\ and p 2 are the initial four-momenta of nucleons 1 and 2, respectively, i.e., the first 
term results from the photon being emitted by nucleon 1 and for the second term, it is being emitted by nucleon 2. 
In the same schematic notation, the four-divergence of the interaction current then is simply [23|, |3Q( 

= VQ-QV 

= Vip'^p'^p! - k,p 2 )Q x + V{p' x ,p' 2 ,pi,p 2 - k)Q 2 - QiV(pi + k,p 2 ;pi,p 2 ) - Q 2 V{p'^p' 2 + k; Pl ,p 2 ) , (20) 

where the arguments of V are nucleon momenta and the momentum dependence of the interaction current is 

V* = V»(k,p , 1 ,p' 2 ;pi,p 2 ) , with p\ +p' 2 + k =p x +p 2 , (21) 

i.e., the momenta pi,p 2 and p[,p 2 are those of the incoming and outgoing nucleons, respectively. Whether the charge 
operators Qi in (|20[) pertain to incoming or outgoing nucleons is clear from where the Qi are placed in the equation. 
In other words, if placed on the right of V, Qi describes the charge of the incoming nucleon i, and if placed on the 
left, it describes the charge of the outgoing nucleon i. Placing the charge operators in this manner is necessary since 
they interact with the isospin dependence of the interaction V. 
The four- divergence of then follows as 

M" = (QG 1 - G^Q) G V + VGO (qGq 1 - G^O) +VQ-QV- VG^QG^ 1 - G^Q)G V 

= -G^QGoV + VGoQGo 1 - VG QV + VQG V . (22) 
For the entire current, we then find 

k^B" = (1 + TG ) (VGoQGq- 1 - G^QGqV - VGqQV + VQG v) (G T + 1) , (23) 
and thus finally, using ([T]), 

fc,5" = TGqQGq 1 - Gq 1 QG q T . (24) 

This is the correct generalized Ward-Takahashi identity (3l| for the bremsstrahlung current providing a conserved 
current for external on-shell nucleons. In a more explicit notation, using the same arguments for T and B^ as for V 
and V^, respectively, in Eqs. (|20|) and (|2T1) . this reads 

k^ = Tljp'^p'^px - k,p 2 )ti(pi ~ k)Q 1 t^ 1 (p 1 ) +T{p' 1 ,p' 2 ;p 1 ,p 2 - k)t 2 (p 2 - k)Q 2 t 2 1 (p 2 ) 

- trVi) Qi Wi + k ) T (Pi + k,p' 2 ;Pi,P2) - t 2 \p' 2 ) Q2 t 2 (p' 2 + k) T{p' x ,p' 2 + k;p x ,p 2 ) . (25) 

I 

The inverse nucleon propagators tj {p) appearing here ensure that this four-divergence vanishes (i.e., that the 
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bremsstrahlung current is conserved) if all external nu- 
cleon legs are on-shell. 



C. Covariant three-dimensional reduction 

To calculate any reaction amplitude in a full four- 
dimensional framework is a daunting numerical task. 
In practical applications of relativistic reaction theories, 
therefore, one often employs three-dimensional reduc- 
tions that eliminate the energy variable from loop inte- 
grations in a covariant manner leaving only integrations 
over the components of three-momenta. Many such re- 
duction schemes can be found in the literature [32H35j |. 
Our results presented below hold true for any reduction 
scheme that puts both nucleons in loops on their respec- 
tive energy shells. 

For hadronic reactions, the primary technical con- 
straint to be satisfied by any three-dimensional reduc- 
tion is the preservation of covariance and (relativistic) 
unitarity. For reactions involving electromagnetic inter- 
actions, there is the additional constraint of gauge in- 
variance. This is a non-trivial constraint since the reduc- 
tion scheme, in general, will destroy gauge invariance as a 
matter of course. Hence, to restore it, one must introduce 
additional current mechanisms as part of the reduction 
prescription for photoprocesses. As we shall see, this can- 
not be done in a unique manner because gauge invariance 
does not constrain transverse current contributions. 

For the NN problem, three-dimensional reductions re- 
sult from replacing the free two-nucleon propagator Go 
by one containing a S function that eliminates the en- 
ergy integration in loops. In the following, we make the 
replacement 



G go 



(26) 



to indicate that the internal integration is a three- 
dimensional one over the three-momentum of the loop. 
To obtain on-the-energy-shell integral equations from the 
Bethe-Salpeter equations ([T]) when using this reduction, 
the external nucleon legs must be taken on shell as well. 
This provides the reduced integral equations 



t = v + vgot = v + tgov , 



(27) 



where lower-case letters v and t (instead of V and T, re- 
spectively) signify that all nucleons — internal and exter- 
nal — are on their energy shell. However, when consid- 
ering below the gauge invariance of the bremsstrahlung 
current that results from the reduction (|2"5)l . we require 
fully off-shell T matrices. They can be obtained from 
iterated versions of (QJ which is then subjected to the 
reduction (|2l)j) . producing 



T = V + V(g + g tg )V 



(28) 



where all external nucleons may be considered off-shell. 
This off-shell T, thus, is obtained by quadratures from 



the integral-equation on-shell solution t. However, this T 
is not the same as the solution of the full four-dimensional 
Bethe-Salpeter scattering equation (jTJ) even though we 
use the same notation to keep matters simple, even at 
the risk of inviting confusion. The T"s appearing in the 
following context always refer to fully off-shell or half off- 
shell versions of the T matrix defined by the quadrature 
formula gSJ. 

The previous proof of gauge invariance of the full 
four-dimensional formalism given in Sec. Ill Bl shows that 
gauge invariance depends on an intricate interplay of 
all hadronic reaction mechanisms. Any approxima- 
tion, therefore, will in general destroy gauge invari- 
ance. Hence, we expect that simply subjecting the full 
bremsstrahlung current of Eq. (|14p to the reduction 
prescription (1261) will not retain gauge invariance, and 
that, therefore, additional steps will be necessary to en- 
sure gauge invariance. Since, from a formal point of 
view, any modification of a given current can, without 
lack of generality, always be expressed by adding an ex- 
tra current, we may write the gauge-invariant current 
that results from a judicious adaptation of the three- 
dimensional reduction procedure to B^ of Eq. (TT4"]) in the 
form 



B^ 1 -> Bl 1 = 



(TG + l)J»(l + G T) 



X^ 



(29) 



] re dj schemati- 



The first term on the right-hand side, 
cally denotes the necessary modifications of the hadronic 
mechanisms in B^ itself and the last term, Xq 1F , is the ad- 
ditional gauge-invariance preserving (GIP) current that 
is to be determined to make B^ gauge invariant. In other 
words, we demand that 



TGqQG 



G 1 QGqT 



(30) 



i.e., that the four-divergence of the reduced current B^ 
remains identical in form to the generalized WTI of 
Eq. (|24l) and we are going to ensure this by choosing 
Xq W accordingly after having determined [■ • • ] rc( j in (|29l) . 
[We repeat here that — completely consistent with the 
covariant three-dimensional reduction — the off-shell T's 
of (f5U|) . and of the subsequent equations in this section, 
are those defined by the off-shell extension ([28]) of ([27]). 
and not the solutions of the original Bcthc-Salpctcr equa- 
tions ([1]) that appear in ([24]) ; see also the discussion sur- 
rounding Eq. pg]).] 

As a first step, we employ the reduction for the external 
Gq factors and write 



Bt 



(Tg + l)J?(l+goT)+X% lp 



(31) 



where is the reduced form of the basic production cur- 
rent of (jlll) . For its determination, we note that one 
cannot simply employ the reduction (]26j) for every Go 
appearing in ([lljl since this produces unphysical mecha- 
nisms. The box-graph contribution 



= VGod^GoV 



(32) 
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V 



V 



V 



V 



V 



(b) 



FIG. 4. (a) Sum of box graphs describing Eq. (J32J, with 
rectangles labeled V subsuming all mechanisms of the NN 
interaction. All internal nucleon lines are off-shell, in general, 
(b) The two possibilities of putting nucleons on-shell in the 
internal loop of the first diagram of (a), as indicated by "x" 
on the corresponding nucleon line, leaving one nucleon at the 
photon vertex off shell. Analogous diagrams can be drawn for 
the corresponding second diagram. 



shown in Fig.SJa), that appears as a subtraction in (fT2]) . 
cannot be reduced in the form 

VG d^G Q V -> Vg (Fg V (unphysical) (33) 

since this would have the bremsstrahlung photon emerg- 



ing from intermediate on-shell nucleons which is not pos- 
sible for a physical photon. At least one of the nucleon 
legs at the photon vertex — either the incoming or the 
outgoing one — must remain off-shell. The two possi- 
bilities of doing that allowing the production of physical 
photons are 



VGodPGoV -> 



f VGo&goV 
{VgodfGoV 



(34) 



as shown in Fig.HJb). Since there is nothing that suggests 
that one choice is to be preferred over the other, we allow 
for both and thus make the replacement 

b» -> K = XiVGod^goV + \fV 9() d»G a V , (35) 

where the constant factors are constrained by 

Ai + A/ = 1 (36) 

to avoid double counting. The symmetric choice would 
be A, = Xf = 1/2, of course, but we want to allow here 
more flexibility for reasons given below. 

In detail, the reduced bremsstrahlung current thus 
reads 



B? = (Tg + 1) [d»G V + VG d» + V" + D?] (1 + g T) + A£ IP , 



where 



Dt l , = n i - W, 



(37) 



(38) 



is the reduced form of (|12|). with denoting the reduced form of the transverse current T M satisfying fc^T/ 1 = 0. 
Note that the Go's appearing in the terms cannot be reduced to go's for the same reason that the reduction (l33|) 
is not possible. Evaluating the four-divergence of this expression provides 

k^B? = (Tg Q + 1) [(QGo 1 - G^Q) G a V + VG (qG 1 - G^Q) +VQ-QV- fe„6Jf)] (1 + g T) + fc M X£ IP 
= -(Tgo + 1)Gq 1 QGqT + TG QG \l + g T) - (Tg + l)k M b?(l + g a T) + k^ lp 



TGqQG q 



G^QGqT - (Tgo + l)fc„6£(l + g T) + k M X^ 



(39) 



It was used here that 

TG o QG o 1 g o T = and Tg GQ 1 QG T — (40) 
vanish identically. For the proof, consider 

goG^QGo -> , Al ^_ 2 , 2 o [t^( P )QMP + k)] , 

s — (ei + £2) 

(41) 

where the right-hand side here provides one generic con- 
tribution (stripped of all extraneous factors) contained 
in the expression on the left. The "o" symbol indicates 
the remaining three-momentum loop integration. The 
variable s is the squared total energy of the system and 



the £j are the individual on-shell energies of the two on- 
shell nucleons in the loop [indicated by the symbol "x" 
in Fig. 2Kb)]. The A; are the positive-energy projectors 
of the nucleons. The energy component of p = (ei,p) is 
on-shell and thus 

A 1 (p)t r 1 ( P )^^-^A 1 (p) = 0. (42) 
Im 

Because of p 2 — m 2 , this term always vanishes for any 
p. Therefore, even if s — (ei + e^) 2 should vanish as well 
for some p, this cannot be compensated, i.e., the limit 
of the corresponding 2 situation is zero. This proves 

that TgoG^QGoT = 0; the proof for TGoQG^gaT = 
follows in a similar fashion. 
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FIG. 5. Reduced basic bremsstrahlung production current of Eq. I|48[l resulting from the covariant three-dimensional 
reduction (see Sec. Ill CP for nucleon-nucleon interactions based on single-meson exchanges. Time proceeds from right to left. 
The mesons M taken into account in the present work and in Ref. [J are tt, rj, p, u, a, and ao (formerly S). The nucleonic 
current corresponds to the first two diagrams on the right-hand side and the meson-exchange current is depicted by the third 
diagram. The first four diagrams respectively labeled s, u, t, and c correspond to the complete gauge-invariant description for 
the process NM — > for the upper nucleon line, with the labels s, u, and t alluding to the kinematic situations described by 
the corresponding Mandelstam variables (see Fig. [6] and the corresponding discussion in Sec. HID) . The fourth diagram contains 
the NM — > N'y four-point contact current M^ t discussed in Sec. Ill El labeled "c" in the diagram. The correct gauge-invariance 
treatment of this diagram was found to be crucial in reproducing the KVI data in Ref. 0]. The diagrams corresponding to s, 
u, and c for the lower nucleon line are suppressed. The last diagram (labeled Jt ) subsumes the transverse five-point currents 
of (|46p . Transverse transition-current contributions are also subsumed in this diagram; examples are the 7p7r, 'yuuir, and jNA 
transition currents depicted in Figs. [3] and [9] respectively. Antisymmetrization of identical nucleons is implied. 



The first two terms in the last line of already pro- 
vide the complete four-divergence (|30p necessary for the 
gauge-invariance condition to hold true. It follows then 
that the last two terms must vanish, 

fc M X£ IP - (T 5o + l)fc^(l + g T) = . (43) 
The four-divergence of Xq IF thus is constrained by 

= (Tgo + i)fc^(i + g T) , (44) 

where describes the longitudinal pieces of the reduced 
box-graph current b% . Hence, without lack of generality, 
we may write 

Xg w + (Tg + l)D?(l+g T) = (Tg +l)J^(l+g T) (45) 

where Jf is the purely transverse current 

J£ = T? - X.VGodt^goV - \ f Vg d%G V ; (46) 

dx here only contains the transverse pieces of the nucleon 
currents as they appear in g? m . As far as gauge invariance 
of Bp is concerned, the current J^ is irrelevant. There- 
fore, the parameters A.; and A/ may now also be treated 
as independent parameters, unconstrained by 

To summarize the present results obtained for the 
three-dimensional reduction, in this approximation the 
gauge-invariant bremsstrahlung current reads 

Bf = (Tg + l)JF(l+g T), (47) 

where the reduced basic production current is given as 

= d"G Q V + VG d» + V + J% . (48) 

How one chooses Jt in an application is not fixed by the 
formalism, beyond the generic form given in (|46[) . Note 
that the generic graphical structure depicted in Fig. [^re- 
mains valid also for J^, with the last graph labeled D 11 on 



the right-hand side of the figure depicting now the trans- 
verse five-point current j£. For the specific case of NN 
interactions based on single-meson exchanges only, the 
reduced basic production current of (T4"5)) is illustrated 
diagrammatically in more detail in Fig. [5] 

Let us add some remarks here. Even though the proce- 
dure to preserve gauge invariance was presented here in 
terms of an additional ad hoc current Xq IP , the derivation 
shows that, rather than adding a current, the application 
of the three-dimensional reduction procedure to the cur- 
rent B^ L really amounts to dropping (at least part of) the 
reduced contribution from 6 M that was necessary in the 
full four-dimensional treatment to prevent double count- 
ing of the TGod^GoT contribution. The particular form 
of of (T4"o| follows from exploiting the constraint (j4"4")l 
to the extent to which it is possible since the gauge in- 
variance cannot constrain transverse contributions. Our 
applications reported here, in Sec. IIII| and in [3] suggest 
that the choice A^ = Xf = yields by far the best numer- 
ical results, in view of the fact that the double-counted 
term TGod^G^T in the full four-dimensional formulation 
of Sec. Ill Al now appears in (|47|) as two distinctly dif- 
ferent contributions TG^d^g^T + Tg d^GoT : there are 
now no longer any terms counted doubly and thus the 
corresponding subtraction is no longer necessary either. 
The finding that A^ — A/ — yields the best numer- 
ical results is completely consistent with this fact. In 
our calculations, therefore, Jt only contains transverse 
contributions from electromagnetic 7p7r, jujit, and 7^A 
transitions (for the latter, see discussion in Sec. Ill F[) . 



D. Local gauge invariance: Constraining 
subprocesses 

The gauge-invariance constraints in terms of the gen- 
eralized Ward-Takahashi identities, either in the form 
for the full amplitude or as (13TJ)) for the reduced 
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one, are just formal constraints, of course, since the only 
physically relevant — i.e., measurable — ramification is 
current conservation when all external nucleons of the 
bremsstrahlung process are on their respective energy 
shells. It is for this reason that the current-conservation 
constraint, 

k^B* 1 = (external nucleons on-shell) , (49) 

called global gauge invariance, is the only constraint 
that is implemented in many reaction models of pho- 
toprocesses. We would like to advocate, however, that 
using this as the sole constraint is not enough if one 
aims at providing a consistent microscopic description 
of the photoreaction at hand. It was shown in the pre- 
ceding sections that the gauge invariance of the total 
bremsstrahlung amplitude hinges in an essential way on 
each subprocess providing its correct current share to en- 
sure the gauge invariance of the entire current — and 
thus ultimately provide a conserved current. This is only 
possible if the current associated with each subprocess 
satisfies its own generalized Ward-Takahashi identity, as 
exemplified here by Eq. (1161) for the propagators and by 
Eq. (|20[) for the interaction current. Thus, imposing local 
gauge invariance, i.e., imposing consistent off-shell con- 
straints of this kind for all subprocesses in a microscopic 
description of the reaction at hand will then automati- 
cally ensure that the total amplitude will satisfy a gen- 
eralized WTI of its own. This not only means that it 
will indeed satisfy the physical constraint of a conserved 
current, but beyond that it will ensure that if the pro- 
cess at hand will be used as a subprocess for another, 
larger reaction, it will automatically provide the correct 
contribution to make the larger process gauge invariant 
as well. 

In this respect, we draw attention to the fact that the 
reaction dynamics depicted by the first four diagrams 
on the right-hand side of Fig. [5] correspond to the cap- 
ture reaction NM —> Nj where the meson M is emitted 
by the spectator nucleon depicted by the lower nucleon 
line. The corresponding time-reversed equivalent meson- 
production process is shown in Fig. [5] for the example of 
pion production. The current amplitude for this process 
can always be broken down into four generic contribu- 
tions (see, for example, [H|), 

= Mjf + M£ + Mf + , (50) 

where the indices s, u, and t allude to the Mandelstam 
variables describing the kinematic situations of the corre- 
sponding diagrams in the figure. The corresponding in- 
teraction current, in particular, is obtained by attaching 
the photon to the inner workings of the meson-nucleon- 
nucleon vertex F according to 

< t = -{iT- (51) 

If we demand now local gauge invariance, the current 
must satisfy an (off-shell) generalized WTI. In addition 




p' P p' P p' P p' P 



FIG. 6. Generic pion photoproduction diagrams 7 + TV — > 
7r + TV. Time proceeds from right to left. If read from left 
to right, this corresponds to the pion-capture reaction as it 
appears in the first four diagrams on the right-hand side of 
Fig.[S]along the upper nucleon line. The last diagram here la- 
beled "c" corresponds to the interaction current arising from 
attaching the photon to the interior of the ttNN vertex. The 
properties of this contact-type four-point current M^ t are es- 
sential to render the entire amplitude gauge invariant, with 
the necessary constraint equation given by (|52[) . At the tree- 
level for undressed hadrons, this term reduces to the Kroll- 
Ruderman term. 

to the trivial contributions resulting from the propaga- 
tor WTIs for the currents associated with the external 
legs of M£, M£, and M t M , similar to (T5JI, this means m 
particular, that the interaction current must satisfy [23J 

fc M M£ t = Q N F U + Q M F t - F S Q N , (52) 

where F x denotes the meson-nucleon-nucleon vertices F, 
with the subscripts x — s,u,t corresponding to Mandel- 
stam variables of the respective kinematical situation of 
the vertices in s-, and f-channel contributions, as de- 
picted in Fig. |H1 Q m and Q n are the charge operators 
of the meson and the nucleon, respectively. Since the 
charge operators interact with the isospin dependence of 
the vertex F, their placement before or after F x is signif- 
icant to ensure charge conservation. Note that Eq. (l52j) 
is the exact analog of Eq. (|2"U)) for the present applica- 
tion, with Eq. (|20[) providing the four- divergence for a 
five-point interaction current for an outgoing photon and 
Eq. (I52p constraining the four-point interaction current 
for an incoming photon. Using the Q notation introduced 
in Eq. (fT7|) . this may be made more obvious by writing 
Eq. ([52) as 

K M L = (Qn + Qm) F - FQ N , (53) 

where the charge operators of the final hadrons appear on 
the left of F and that of the initial hadron on its right, as 
it is appropriate for the interaction current of any meson- 
production process 7TV — > TVM similar to Fig. H>] for 
the reverse process TVM — > TV7, one needs to change 
k —k and exploit the isospin dependence of the vertex 
to write QmF — —FQm since the meson changes from 
being outgoing to incoming. 

Demanding consistency of the microscopic dynamics 
across various reactions means that the entire TVM — > 
TV7 subprocess in Fig. [5] must satisfy the same general- 
ized Ward-Takahashi identity as the amplitude M M itself, 
except for trivial modifications arising from the fact that 
this is the time-reversed process. This must be true for 
any one of the exchanged mesons — whether scalar, pseu- 
doscalar, or vector — in the bremsstrahlung process, not 
just for the pion example shown in Fig. [5] 
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FIG. 7. Generic structure of interaction current V£sc, Eq. (f55)l . for the single-meson exchange contribution Vuec of the iViV 
interaction depicted on the left. Time proceeds from right to left. The three diagrams on the right-hand side comprise the 
meson-exchange current and the two interaction-current contributions M?, i = 1,2, where the photon interacts with the 
interiors of the respective vertices for nucleons 1 and 2. 



To illustrate in more detail how the requirement of 
local gauge invariance ties together the various current 
mechanisms, let us consider a single-meson-exchange con- 
tribution Vmbc to the full NN potential V. Generically, 
we may write 

Kiec = F 1 t M F 2 , (54) 

where the Fi are the meson-nucleon-nucleon vertices for 
nucleon £ = 1,2 (including all coupling operators and 
isospin dependencies) and tu describes the propagator 
for the exchanged meson M. Graphically, the process is 
given in the diagram on the left-hand side of Fig. If we 
now attach an outgoing photon to all internal mechanism 
of V MEC , we obtain the three diagrams on the right-hand 
side of the figure that comprise the corresponding AW 
interaction current, V^ EC — — {Vuec} 11 ■ This current may 
be written as 

V&a = -F x {LmY F 2 - {Fx}" t M F 2 - Fxt M {F 2 }" 
= F^M T M t M F 2 + M?t M F 2 + Fxt M M% , (55) 

where Mf* — — {Fi}" is the four-point interaction current 
for the vertex £ = 1,2 and 

-{t M Y = t M v M t M (56) 

produces the current operator for the exchanged me- 
son that satisfies the single-particle WTI [28, 29] 

k » T M = <9m*m - ^jQm ■ (57) 

Now, assuming without lack of generality that four- 
momentum and charge flow from vertex 2 to vertex 1 
in Fig. [71 the analogs of (1551) read 

fc^Mf = Fx (Qx + Qm) - QiFx , (58a) 

^M 2 " = F 2 Q 2 - (Q M + g 2 ) F 2 , (58b) 

and thus 

^m^mec — FxtM QmF 2 — FxQMtuF 2 

+ F x (Qx + Qm) t M F 2 - QxFxt M F 2 
+ Fxt M F 2 Q 2 - Fxt M (Qm + Q 2 ) F 2 

= F\tMF 2 Qx — QxFxtMF 2 

+ FxtMF 2 Q 2 — Q 2 FxtMF 2 
= V MEC Q - QV MEC (59) 



which is precisely the gauge-invariance constraint (I20[) 
for this particular contribution to F. In this consis- 
tent microscopic treatment of all subprocesses, therefore, 
the constraints local gauge invariance places on the four- 
point interaction currents Mft t of meson-production pro- 
cesses translate seamlessly into the corresponding con- 
straints on five-point currents of the bremsstrahlung pro- 
cess. The essential step here is to ensure the validity of 
(|58ap for the four-point current in the diagram labeled 
"c" in Fig. [5J and of (j58b[) for its counterpart for the 
lower nucleon line (not shown in Fig. [S]). These relations 
must be true for each of the exchanged mesons. 



The problem of how to ensure the validity of the con- 
straints (JS^J) or (1551) has been studied extensively by the 
present authors and their collaborators for the equiva- 
lent photoproduction processes, and based on the origi- 
nal ideas presented in Refs. [23, 25] , a general prescription 
was given in Ref. [24| that is applicable just as well for 
the most general case of explicit final-state interactions 
with completely dressed hadrons as it is for phenomeno- 
logical vertex functions. The necessary extension in the 
context of bremsstrahlung — to account for the virtual 
nature of the incoming and outgoing nucleons and the 
exchanged meson at the four-point vertex — is accom- 
plished following the work of Ref. [36| . 



For our present applications, and those of Ref. [|[, 
based on single-meson exchanges containing phenomeno- 
logical meson-nucleon-nucleon form factors, the details of 
the four-point currents Mf and M% and of the complete 
descriptions of the corresponding diagrams depicted in 
Fig. UJ were already given in Ref. J3J; we will not repeat 
them here. We emphasize, however, that the inclusion 
of this interaction- type current, with the correct gauge- 
invariance dynamics that ensure the validity of Eqs. (|58|). 
is essential to bringing about the good quality of our 
bremsstrahlung results, in particular, the good descrip- 
tion of the high-precision KVI data [l|, [J reported in 
Ref. Q. This is a novel feature of our approach that re- 
solves a longstanding discrepancy between the data and 
their theoretical description. 
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E. Interaction current V** for phenomenological 
NN interactions: A systematic gauge-invariance- 
preserving approximation 



The previous discussions show that the gauge- 
invariance conditions (fSTTj) and (1551) for the respective in- 
teraction currents are of crucial importance to ensure an 
overall bremsstrahlung amplitude that satisfies gauge in- 
variance. The understanding here is that the detailed re- 
action mechanisms that go into providing the details of 
these interaction currents are completely known and that 
therefore it would be straightforward to make sure that 
the corresponding gauge-invariance conditions are indeed 
satisfied. This is generally the case only for NN interac- 
tions based on single-meson exchanges between nucleons 
without any phenomenological form factors. And if one 
employs phenomenological form factors, one must resort 
to the prescriptions given in Refs. [2314251 [36| to ensure 
Eqs. (f58|) . as discussed in the preceding section. How- 
ever, beyond problems associated with phenomenological 
form factors, it is conceivable that some details would not 
be available in some instances either because the micro- 
scopic coupling of the photon to the internal dynamics 
of the interaction would not be feasible or sensible or 
because it would be too complicated for practical appli- 
cations. Examples of the former would be NN interac- 
tions based on position-space methods where there are no 
(momentum-dependent) exchange mechanisms that per- 
mit the coupling of a photon in a dynamically meaning- 
ful way. Examples of the latter might be baryon con- 
tributions beyond the nucleon since they require box- 
type iViV contributions with intermediate non-nucleonic 
baryonic contributions that might be too cumbersome to 
be treated with explicit photon couplings (see Sec. Ill Fl 
for A A and/or A A systems and the resulting graphical 



structures depicted in Fig. [8]). 

Nevertheless, each of the corresponding interaction- 
current contributions must satisfy its appropriate gauge- 
invariance condition otherwise the gauge invariance of the 
entire amplitude breaks down. We will show here how 
one can ensure that gauge invariance is preserved for any 
phenomenological contribution to the NN interaction by 
constructing a phenomenological interaction current that 
satisfies appropriate constraints. 

To this end let us assume the total NN interaction V 
can be split up into n independent contributions, 

V = Vi + V 2 + ■ ■ ■ + V n . (60) 

Formally coupling a photon to each contribution accord- 
ing to = — {Vi}, the total interaction current then 
breaks down accordingly into n independent contribu- 
tions, 

v = vj* + y 2 " + • ■ ■ + vj? (6i) 

which, because of their independence, must each sepa- 
rately satisfy a gauge-invariance condition similar to (|20p , 
i.e., 

h^Vf = ViQ-QVi , i = l,2,...,n, (62) 

otherwise the total current could not satisfy lj2"0|) . 

Let us assume now that one of the Vi is such that the 
construction V? = — {V}^ is not readily available. In 
this case, we may devise an auxiliary phenomenological 
current Vf instead that satisfies exactly the same four- 
divergence relation (|B"2l as it would have to be satisfied 
by V^ if it were available. To this end, we adapt the 
procedure used successfully in Refs. [23j-|25| for the four- 
point interaction currents Mf^ of meson production to 
the present five-point-current case and make the ansatz 



VI 



Vf 



V l (p' l ,p' 2 \Pi - k,p 2 ) - Wi(k,p' 1} p2;pi,p2) 



V i (p' 1 ,p' 2 ;pi,P2 -k)- W i {k,p' 1 ,p' 2 \pi,p 2 ) 
Qi(2p£ + A0" 



Qi(2 P i-fc)" 
(pi - k) 2 - p\ 

Q 2 (2 P2 - ky 



(p' 1+ k) 2 -p' 2 

Q 2 {2 P ' 2 + ky 



(pa - k) 2 - p\ 
Vi(p' x + k,p' 2 ;p!,p 2 ) - W i (k,p' 1 ,P2',Pi,P2) 



(p' 2 + k) 2 -p' 2 
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Vi(p[,p' 2 + k;pup 2 ) - W l (k,p' 1 ,p 2 ;p 1 ,p 2 ) 



(63) 



where W% is a function to be chosen to ensure that each term here is free of propagator singularities. The four- 
divergence of this auxiliary current is then readily seen to produce indeed 



KVt = VtQ- QVi 



(64) 



W i (k,p' 1 ,p' 2 ;pi,p 2 ) (Qi + <92) init iai - (Qi+Q2) Rna iW i (k,p' 1 ,p' 2 ; Pl ,p 2 ) = 



(65) 
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vanishes because of charge conservation. The ansatz (l63l) thus preserves the gauge-invariance condition for the entire 
interaction. As to how to choose Wj, one of the simplest possibilities is 

W i (k,p' 1 ,p' 2 ;pi,P2) = V l {p' 1 ,p' 2 ;pi,p 2 ) 

+ [{p l -k) 2 -p\] [{ P2 -kf-pl] [(pi + fc) 2 -Pi 2 ] [(pi + kf-p'^Riik^p'^pz) , 

(66) 

where Ri(k, pi, p' 2 \ p\, p 2 ) , except for symmetry constraints, is largely arbitrary (and may be equal to zero). This 
choice ensures that the resulting GIP current is free of kinematic singularities, as it must be. In detail, one thus has 



T/ /x Vi(p'l,P2>Pl ~ fc >P2) ~ Vi(p' 1 ,P2,Pl,P2) „ ,„ ,sa Vi(p' 1 ,P 2 ,Pl,P2 ~ k) - Vj (pj , p' 2 j gi ,p 2 ) n , 

1 { Pl -kf-p\ Qi(2P!-fc) (pz-kf-pl Q2(2p 2 -k) 

n fc, f , w Y*M ± k 'P2iPl:P2) - Vi(p' 1 ,p' 2 ;Pl,P2) n , , ,„^ ^(Pl,P2 +fc;Pl,P2) - Vj(pj,P2;Pl»P2) 

- Ql(2Pl + fc) (Pi + A) 2 -Pi 2 P2 + } (p 2 + ^-p 2 2 

+ i?,(fc,pi,p 2 ;p 1 ,p 2 ) [(pi + fc) 2 - pi 2 ] [(p' 2 + fc) 2 -p' 2 2 ] 

x { [(pa - fc) 2 -p 2 ] Q!(2p! - ky + [(px - fc) 2 -pi] Q 2 (2p 2 - fc^} 

+ {g!(2pi + ft)" [(p' 2 - fc) 2 - p' 2 2 ] + g 2 (2p 2 + kf [(Pi - kf - pi 2 ] } 

x [( Pl - k) 2 - pj] [(p 2 - fc) 2 - pi] i? 4 (fc,pi,p 2 ;p!,p 2 ) . (67) 



Note that the subtracted potential contribution 
Vi(pi,P 2 ;pi,P2) is unphysical since pi +p 2 ^ pi +p 2 . 

Hence, the procedure just outlined allows a systematic 
hybrid treatment of all independent contributions Vi to 
the full interaction V where some five-point currents V? 
can be treated explicitly and some may be replaced by 
auxiliary currents Vf according to (|63|) . The total in- 
teraction current V 1 will satisfy the condition (|2"0)l as a 
matter of course and gauge invariance is not at issue. 



F. Extension to coupled channels: NN, NA., and 

AA 

The derivation of the bremsstrahlung current in 
Sec. Ill Al is completely generic and will remain true re- 
gardless of the actual mechanisms taken into account in 
the nucleon-nucleon interaction V that drives the Bethe- 
Salpeter equation ([1]). The current mechanisms depicted 
in Fig. [3] assume an interaction based on single-meson 
exchanges between nucleons only. Here we briefly dis- 
cuss the necessary modifications if such exchanges involve 
transitions between different baryonic states. We limit 
the discussion to transitions between the nucleon A and 
the A(1232) mediated by single-meson exchanges, i.e., 
we consider the effect of coupling the channels AA, NA, 
and AA; transitions into other (resonant) baryonic states 
can be treated along the same lines. 

It is a very simple and straightforward exercise to de- 
couple the corresponding set of Bethe-Salpeter equations 



that couple the NN, NA, and A A channels and write 
the NN interaction appropriate for the smg/e-channcl 
Bethe-Salpeter equation (JlJ as 

V = V UEC + V A ■ (68) 

The first term, V MEC , describes single-meson exchanges 
between nucleons, as shown on the left-hand side of 
Fig. that provide the current mechanisms depicted on 
the right-hand side of the figure. The second term, Va, 
contains all intermediate AA or AA contributions and 
their transitions. Using the notations 

Un- meson-exchange transition A A — ► A A (69a) 

and 

Ua- meson-exchange transition A A — > A A (69b) 

for the transition interactions that mediate the coupling 
to the primary AA channel, we obtain 

Va = (Gpj + G^T^ N G^ L/jv + U A G A T AN G^U jv 

+ U A (G A + G A T AA G A ) Ua + U n G^T0 a G a Ua 

(70) 

where G^ and G A describe the intermediate propagation 
of the AA and the AA systems, respectively. Interme- 
diate transitions AA -> AA, AA -> AA, AA -t AN 
AA — > AA are subsumed in the respective T matrices 

rpA rpA rpA „ J T^A 

1 N,N^ A,JV> AT, A' dllu J A,A- 




FIG. 8. Contributions to the NN interaction involving inter- 
mediate iVA and AA systems and intermediate transitions 
Tj\ between such systems as given by Va in Eq. ((70)) . Sum- 
mations over all possible exchanged mesons that mediate tran- 
sitions to the A state are implied. Here, we take into account 
7r and p exchanges. Antisymmetrization of external nucleons 
is implied. 

Figure [5] provides a graphical representation of Va- 
Coupling a photon to each of the mechanisms depicted 
here results in the interaction current 

vi = -{v A r (71) 

which, together with the interaction current Vsiec de- 
picted generically in Fig. [Jj constitutes the total inter- 
action current 

V = V&c + VI (72) 

(if nucleons and A's are the only baryon degrees of free- 
dom considered). Of course, calculating V£ explicitly is 
a formidable task. There are eight current contributions 
for each of the simple box graphs (where the photon can 
couple to each of the four intermediate hadrons and the 
four vertices) and eleven for each of the graphs involv- 
ing an intermediate T matrix. Considering only tt and p 
exchanges for these graphs involving the A and, allow- 
ing for the four possibilities of exchanging these mesons, 
altogether, therefore, there are 240 current terms. This 
is without separately accounting for the internal mecha- 
nisms resulting from coupling the photon to the T ma- 
trices with A degrees of freedom which is similar in com- 
plexity to the NN bremsstrahlung current itself. In view 
of this complexity, we forego drawing the corresponding 
diagrams for V£. 

Each set of currents resulting from one graph in Fig. [8J 
corresponds to an independent interaction current, as dis- 
cussed in conjunction with (|6ip . and therefore must sat- 
isfy the gauge-invariance condition (|62[) independently. If 
an exact treatment of the corresponding interaction cur- 
rent is not feasible in view of the complexity of the prob- 
lem, from the point of view of gauge invariance this may 
also be done in an approximate manner by constructing 
an auxiliary current for each of the graphs in Fig. [8] along 
the lines discussed in the preceding Sec. Ill El 

In addition to the hadronic transitions into intermedi- 
ate A states, there are also direct electromagnetic tran- 
sitions 7./V — > A. Such contributions cannot be ob- 
tained by applying the gauge-derivative method used in 
Sec. Ill Al or by any other method based on minimal sub- 
stitution. They must be added by hand, in terms of their 
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FIG. 9. Lowest-order electromagnetic 7iVA-transition contri- 
bution to of Eq. (|73[) ; the exchanged mesons are M — tt, p. 
The full contribution involves replacing the single-meson ex- 
changes in the two graphs by the T-matrix representing the 
sum of all two-nucleon irreducible transitions NN —¥ AN and 
AN — > NN, respectively. The transition currents are trans- 
verse and thus have no bearing on gauge invariance. Anti- 
symmetrization of nucleons is implied. 

own Lagrangian [see Appendix, Eq. (|A.lc|) ]. Their full 
contribution to the basic production current J p of (fTT|) 
may be written as 

J a = dpjA Gn Tan,nn + Tmn,an d AN , (73) 

where Tan, aw and Tnn,an are the T- matrices result- 
ing from summing all two-nucleon irreducible transitions 
NN ->• AN and AN -> NN, respectively; d^ A and 
d AN , in an obvious schematic notation borrowed from the 
NN contributions, contains the electromagnetic transi- 
tion current A — ► N and N — > A, respectively, along 
the baryon lines that emit the photon, as depicted in 
Fig. [9] for the lowest-order single-meson exchange con- 
tributions to the respective T-matrices. These contribu- 
tions are manifestly transverse and hence have no impact 
on gauge invariance. Graphically, they may be subsumed 
in the right-most diagrams of Figs. [2] or [5j depending 
on whether one considers the full formalism or its three- 
dimensional reduction, respectively. In other words, 
given above is part of the transverse current T M , as an- 
ticipated already in Eq. (p~5)) . 

In summary, the full basic interaction current for AW 
bremsstrahlung, including A degrees of freedom, be- 
comes 

J" = d^GoV + VG Q d» + V» + [J& + J£ - VG d^G V] , 

where the groupings of the four terms corresponds to the 
four graphs on the right-hand side of Fig.[2j Here, V and 
V* contain A degrees of freedom according to Eqs. (155)) 
and (|72l) , respectively, and the transverse current is taken 
as T M = J^ f + J^. In the reduced case, we have 

J? = #G V + VG # + V' 1 + [J^ + J£J , (75) 

where the last grouping is the explicit transverse current 
Jt = J^j + J^a> w ith Ai and A/ of Eq. (|4"o]l put to zero; 
the reduced N A-transition current J^ A is obtained by the 
corresponding three-dimensional reductions of the loop 
integrations within J^. 

In our present results (discussed in the next section), 
in view of their complexity, we do not take into account 
any currents that account for the box-graph mechanisms 
of Fig. [SJ i.e., V£ is set to zero. However, we do take 
into account the transition-current contributions (|73|) in 
lowest order, as depicted in Fig. |H1 
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III. APPLICATIONS 

The present formalism was applied successfully Q in 
the description of the high-precision KVI data [lj at a 
proton incident energy of 190 MeV. In this case, in view 
of the relatively low energy, the dominant current con- 
tribution is expected to be essentially nucleonic. We 
felt justified, therefore, to restrict ourselves to the NN 
channel only and neglect the effects of coupling to other 
baryon-baryon channels. Here, we apply the formalism to 
some selected data sets from TRIUMF [26[ at a higher 
incident energy of 280 MeV. This energy is just above 
the pion-production threshold energy and still well be- 
low the A(1232) resonance peak energy region of about 
650 MeV. Nevertheless, an earlier analysis of Ref. [l8[ 
within a coupled-channel approach (employing NN and 
NA channels) indicated that incorporating A contribu- 
tions leads to better agreement with the data at certain 
geometries. For the TRIUMF experiment at 280 MeV, in 
particular, the dominant A contribution was found to be 
arising from the jNA transition current, J^, in lowest 
order [18J (a similar result was found in the single-channel 
approach of Ref. [Hj|). We point out, however, that the 
approach reported in Ref. [181 ] is not gauge invariant be- 
yond the soft photon approximation. Since here we wish 
to preserve full gauge invariance, we need to follow the 
procedure outlined in Sec. Ill Fl for the incorporation of 
A degrees of freedom. In general, however, this is too 
challenging technically at present, except for the lowest 
order of the transition-current contributions of Eq. (|T3[) 
depicted in Fig. [9J which we do take into account in the 
results presented below. Incorporation of the interaction 
current Vj? arising from the A-box-diagram contributions 
of Fig. [8] is beyond the scope of the present paper and we 
will leave this for future studies. 

The details of the present model for constructing the 
production current J M as well as the NN interaction are 
given in the Appendix. 

Figure [TU] shows a comparison of the present results for 
the cross sections with the TRIUMF data in the coplanar 
geometry at small proton scattering angles of 9\ = 12° 
and 6*2 = 14° and at larger proton scattering angles of 
B\ = 28° and 62 = 27.8°. We mention that in the present 
work, we follow the convention of Ref. [3, Q for assigning 
the nucleon angles 9\ and 9 2 which is just opposite to 
the one adopted in other studies [H, Q11 [Hj]. Also, no 
arbitrary normalization factor of 2/3 to the data points 
[2(^| is applied here. The parameters h and a of the phe- 
nomcnological four-point contact current determined in 
Ref. [3] in the description of the KVI data are kept at 
the same values, h = 2.5 and a — 1000, despite the fact 
that they had been determined at the much lower inci- 
dent energy of 190 MeV. In principle, of course, these 
parameters may be functions of the incident energy. We 
see that the agreement is very reasonable. The influence 
of the A-resonance contribution is relatively small but 
significant and definitely helps to improve the agreement 
with the data. This is in line with the earlier findings 
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FIG. 10. (Color online) Comparison of the present prediction 
for the cross sections with the TRIUMF data [26} in coplanar 
geometry at 280-MeV proton incident energy. 81 and 82 de- 
note the fixed scattering angles of the two protons in the final 
state. The cross sections are shown as functions of the emitted 
photon angle, # 7 , in the laboratory frame. The (blue) solid 
curves show the results with the total current of Eq. (|75|) . 
The (red) dashed curves correspond to the results when the 
transition current, = J^j + J£a> ^ s switched off in its en- 
tirety. The transition-current contributions are shown sepa- 
rately as (green) dash-dotted curves (which are hardly visible 
in the graphs) for the meson transitions (771 p and -ynco) sub- 
sumed in and as (brown) dot-double-dashed curves for 
the 7 A A transitions that provide the lowest-order contribu- 
tions of J^ A . The (magenta) dotted curves correspond to the 
results when the generalized four-point contact current (cf. 
Fig. [5]) is switched off. 



of Refs. [ljl fl8| . Furthermore, the present results are 
of comparable quality to the earlier results [r| [l8[ in 
reproducing the TRIUMF data, indicating that, unlike 
for the highqjrecision KVI data at small proton scatter- 
ing ang les [JH, the generalized four-point contact cur- 
rent — required to maintain the gauge invariance of the 
full bremsstrahlung amplitude — is apparently not that 
critical. The dotted curves in Fig. [TU] correspond to the 
results when the generalized four-point contact current 
is switched off. As one can see, its effect is visible but 
somewhat smaller than the effect of the A contribution. 
Overall, however, apart from the differences in reaction 
energies and scattering-angle configurations, we note that 
the TRIUMF data are much less precise than the KVI 
data, thus making the assessment of how relevant vari- 
ous theoretical contributions are in reproducing the data 
much less conclusive. 
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FIG. 11. (Color online) Comparison of the present prediction 
for the analyzing power with the TRIUMF data [2(| in copla- 
nar geometry at 280-MeV proton incident energy. 9i and 62 
denote the fixed scattering angles of the two protons in the 
final state. The analyzing power is shown as a function of the 
emitted photon angle, # 7 , in the laboratory frame. The (blue) 
solid line is obtained with the total current J^ 1 of Eq. (|75|) . 
The (red) dashed curve corresponds to the result when the 
7./VA transition current, J^ A , is switched off. 

In Fig. 1 1 1 1 we show our results for the analyzing 
power together with the TRIUMF data at 6 1 = 12° and 
82 = 12.4°. The agreement is reasonable given the un- 
certainties in the data. Again, the A-resonance current 
helps improve the agreement with the data. 

In summary, therefore, our results for the TRIUMF 
data at a proton incident energy of 280 MeV are in agree- 
ment with the findings of the previous studies [l6llT8|. 
We point out, however, that our results are obtained 
within a fully gauge-invariant framework. Moreover, one 
may hope that inclusion of the A-box contributions may 
further improve the agreement. 



IV. SUMMARY 

We have presented a complete, rigorous formulation 
of the NN bremsstrahlung reaction based on a relativis- 
tic field-theory approach in which the photon is coupled 
in all possible ways to the underlying two-nucleon T- 
matrix obtained from the corresponding covariant Bethe- 
Salpeter-type NN scattering equation using the gauge- 
derivative procedure of Haberzettl J23J. The result- 
ing bremsstrahlung amplitude is unitary as a matter of 
course and it satisfies full local gauge invariance as dic- 
tated by the generalized Ward-Takahashi identity. The 
novel feature of this approach is the consistent — i.e., 
gauge-invariant — incorporation of interaction currents 
resulting from the photon coupling internally to interact- 
ing hadronic systems. 

The formalism is quite readily adapted to approxima- 
tions and thus can be applied even in cases where the 



microscopic dynamical structure of the underlying inter- 
acting hadronic systems is either not known in detail or 
too complex to be treated in detail. We have pointed 
out how the interaction currents resulting from the pho- 
ton being attached to nucleon-nucleon-meson vertices can 
be treated by phenomenological four-point contact cur- 
rents that preserve gauge invariance following the ap- 
proach of Haberzettl, Nakayama, and Krewald [2J|. In 
an advance application of the present formalism [3[ , such 
interaction currents had been shown to contribute signif- 
icantly to reproducing the high-precision proton-proton 
bremsstrahlung data at 190 MeV obtained at KVI |l), 
thus removing a longstanding discrepancy between the- 
ory and experiment. In addition, we have provided a 
scheme that permits the approximate treatment of cur- 
rent contributions resulting from pieces of the NN inter- 
action that cannot be incorporated exactly. In each case, 
the approximation procedure ensures gauge invariance of 
the entire bremsstrahlung amplitude. 

We have also discussed the necessary modifications 
when taking into account baryonic states other than the 
nucleon N; in detail, we consider the A(1232) resonance 
by incorporating the couplings of the NN to the N A and 
AA systems, and the 7 — > NA transitions. 

The formalism has been applied to the 280-MeV 
bremsstrahlung data from TRIUMF (2(| and we find 
good agreement with the data similar to what was ob- 
tained by other authors. We emphasize, however, that 
our results obey full gauge invariance, whereas previous 
approaches provide a conserved current only in the soft- 
photon approximation. Since the gauge-invariant incor- 
poration of intermediate NA and AA configurations is 
too demanding technically at present, the only A de- 
grees of freedom employed in this calculation have been 
the 7./V A-transition currents. 

Finally, we emphasize that despite its completeness 
and generality, the present approach is quite flexible and 
amenable to approximations, as discussed in Sees. Ill Dl 
and III El We expect, therefore, that it will also be use- 
ful in the on-going investigation of hard bremsstrahlung 
d, Q, as well as in di-lepton production processes 0| 
where the photon is virtual. 



Appendix 

For the present application, we construct the basic pro- 
duction current in the framework of a covariant three- 
dimensional Blankenbecler- Sugar reduction [32j (see also 
Ref. [33j]), as described in Ref. H, consistent with the 
NN FSI and ISI based on the OBEP-B version of the 
Bonn potential [37] . In addition to the resulting produc- 
tion current used already in Ref. Q , the present work 
includes the A-resonance current interaction calculated 
in terms of the following Lagrangians: 

Can, = %^A^f t . (d^)N + H. c. , (A.la) 



1G 



C-ANp 



Canj 



(2) 

AM 2 '5 P^ 
(3) 

9an 



AMI 



f A^gf t • (d v p„ v ) N + H. c. , (A.lb) 



o (1) 



2M 

(2) 



(A.lc) 



A7 



= 2.12 

(3) 



where = c^p,, - d v p^ and i* 1 ^ 
The coupling constants are taken to be <7ajv 

an d 9a Np = 40 - 5 Si; for simplicity, g£> Np and g£ Np 
are set to zero. The electromagnetic coupling constants 

9AN-y — 5-00 an d 9an-i = — 4.73 are extracted from the 
corresponding helicity amplitudes [Hj]. The NNtt and 
NNp vertices are derived from the corresponding La- 
grangians used for constructing the nucleonic current in 



Ref. 0. 

Also for simplicity, the AY7T and ANp vertices are 
provided with the same form factors as in the NNtt and 
NNp vertices Q for the off-shell meson and baryon. 

For the A propagator, we use 



i> + Ma 



9 2 -M A +i^ 



3M A 



+ ;rW + 



3M 



Li V 

2~P P 



where M A = 1232 MeV and T A = 120 MeV. 



(A.2) 
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